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Substantial amounts of the transition metals Mn, Fe, Co, and Ni can be substituted for Li in 
single crystalline Li2(Lii_2, r2,)N. Isothermal and temperature-dependent magnetization measure¬ 
ments reveal local magnetic moments with magnitudes significantly exceeding the spin-only value. 

The additional contributions stem from unquenched orbital moments that lead to rare-earth-like 
behavior of the magnetic properties. Accordingly, extremely large magnetic anisotropies have been 
found. Most notably, the magnetic anisotropy alternates as easy-plane —^ easy-axis —>■ easy-plane —>■ 
easy-axis when progressing from T = Mn —> Fe —^ Co —>■ Ni. This behavior can be understood based 
on a perturbation approach in an analytical, single-ion model. The calculated magnetic anisotropies 
show a surprisingly good agreement with the experiment and capture the basic features observed 
for the different transition metals. 


Magnetic anisotropy or magnetic anisotropy energy 
(MAE) is a fundamental concept in solid state science 
affecting magnetic data storage, permanent magnets and 
the investigation of various basic model systems. In a 
simple picture, MAE is the energy necessary or reori¬ 
ent magnetic moments in a certain material. It’s value 
is largely determined by the single-ion anisotropy of the 
magnetic centers. This anisotropy stems from the or¬ 
bital contribution to the magnetic moment (either di¬ 
rectly or indirectly via spin-orbit coupling). Significant 
orbital contributions to the magnetic moment and the re¬ 
sulting large MAEs are usually associated with rare-earth 
elements. In contrast, the orbital moment in 3d transi¬ 
tion metals is normally quenched by the crystal electric 
field. Accordingly, the magnetic anisotropy of these ele¬ 
ments is often small or non-existent. 

Recently, we have found a remarkable exception to 
this trend: iron, when substituted in lithium nitride, 
Li 2 (Lii_a;Fea;)^ behaves in many aspects like a rare- 
earth element^. With an estimated MAE of several 
hundred Kelvin and, in accordance, an observed coer- 
civity field of more than 11 Tesla this compound ex¬ 
ceeds even the largest values observed in rare-earth based 
permanent magnets. Besides iron, other 3d transition 
metal substitutions were synthesized, in polycrystalline 
form, as early as 1949 by Sachsze and duza^^ and have 
been subjects of ongoing experimental^^^^ and theoret¬ 
ical investigations^iffiS. It has been found that T = 
Mn, Fe, Co, Ni, and Cu can be substituted for one of the 
Li-sites: the two-fold coordinated lb Wyckoff site. In¬ 
dications fora^ unusual oxidation state of -1-1 were also 
reportedl^lSHZliil transition metals carry a sizable 

local magnetic mom ent e xcept for T = Cu and highly 
concentrated T — Due to a lack of large enough 

single crystals, there has been no direct access to the 
anisotropy of the physical properties. Only recently we 
developed a single crystal growth technique that is based 
on a lithium-rich flux and is applicable to Li 2 (Lii_ 2 : Ta,)N 
as well as other nitrides and lithium based compounds^ 

Here we present the magnetic anisotropy of 
Li 2 (Lii_a;T 2 ,)N for T = Mn, Co, and Ni and com¬ 


pare the results to our earlier T = Fe work. Two 
concentrations, a dilute one and a more concentrated 
one, were grown, under similar conditions, for the differ¬ 
ent transition metals. The starting materials were mixed 
in a molar ratio of LiiTiLiaN = 8.97:0.03:1 and 8.7:0.3:1 
for dilute and more concentrated samples, respect ively. 
The mixtures were packed in a three-cap Ta crucible^^^^^ 
heated from room temperature to T = 900° C over 5h, 
cooled to T = 750°C within 1.5 h, slowly cooled to 
T = 500°C over 60 h, and finally decanted to separate 
the single crystals from the excess flux. A picture of 
three representative single crystals on a millimeter grid 
is shown in Fig. below. The actual transition metal 
concentration (as opposed to the nominally melt com¬ 
positions given above) in the investigated single crystals 
was determined by chemical analysis using inductively 
coupled plasma mass spectrometry (ICP-MS)®. The 
ICP-MS instrument was provided by Analytik Jena. The 
deviation from the initial concentration, with respect to 
nitrogen, differs from one transition metal to the other 
and also depends on the concentration (see Fig.[^ for 
the measured transition metal concentrations). How¬ 
ever, the obtained concentrations, x, clearly reflect the 
different initial values and allow us to study the dilute 
and concentrated regimes. Magnetization measurements 
were performed using a Quantum Design Magnetic 
Property Measurement System equipped with a 7 Tesla 
magnet. The MAE was calculated analytically in a 
single-ion model based on second order perturbation 
theory using the Green’s function method. 

Figure[2 shows the isothermal magnetization in Bohr 
magnetons per transition metal ion at T = 2 K. The mea¬ 
surements were performed for magnetic field applied par¬ 
allel and perpendicular to the hexagonal c-axis, shown 
by black squares and red circles, respectively. Results 
obtained for the dilute transition metal substitutions are 
shown on the left-hand side, more concentrated substitu¬ 
tions are shown to the right. The largest available field 
of tqH = 7 T allows for saturation of the magnetization 
only for Fe and dilute Ni (and comes close to saturation 
for dilute Mn). Furthermore, the anisotropy field (cross- 
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FIG. 1. (color online) Isothermal magnetization loops at 
T = 2 K in Bohr magnetons per transition metal ion (FC 
±7T refers to a field cooled measurement in +7 T and -T, 
respectively). The sign of the magnetic anisotropy changes in 
an alternating fashion from easy plane to easy axis indepen¬ 
dent from the concentration of the transition metal. 


ing point of the M-curves for H \\ c and _L c) lies well 
above the largest available fields. Therefore, it is not 
possible to accurately quantify the MAE from our data 
(for T 7 ^ Fe). However, the alternating change from easy- 
plane to easy-axis behavior is evident and independent of 
X. Further trends can be recognized: (I) The anisotropy 
decreases with increasing x, except for T — Fe. (II) The 
magnetization values that are approached for /Xpi? = 7 T 
decrease with increasing x. (Ill) The anisotropy observed 
for T = Mn is significantly smaller than that observed 
for the other transition metals (except for concentrated 
Ni). (IV) Even though the sign and magnitude of the 
anisotropy of dilute Fe and Ni appear to be very similar, 
the large hysteresis found for Fe-substitution is absent 
for Ni. (It is also absent for the planar T = Mn and Co.) 
Demagnetization fields amount to only a small fraction 
of the applied magnetic fields and can be neglected. 

Further information about the orbital moment con¬ 
tributions can be obtained from the effective magnetic 
moments, determined from the temperature-dependence 
of the magnetization. Figurej^ shows the temperature- 
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FIG. 2. (color online) Inverse magnetic susceptibility per mol 
transition metal as a function of temperature (x~^ = H/M, 
fj,oH — 7T). A pronounced magnetic anisotropy is observed 
over the whole temperature range except for T = Mn and 
concentrated, T = Ni Li 2 (Lii_a; T' 2 ,)N. The solid lines given 
for the dilute case (left panel) show the inverse susceptibility 
after subtracting the core diagmagnetism of LisN. 


dependent magnetic susceptibility, x(T) = M/H, for 
both orientations of the applied field. Reasonable agree¬ 
ment with Curie-Weiss behavior is observed for T > 
150 K. The contribution of the core diamagnetism of the 
LisN host significantly affects the diluted systems but can 
be neglected for the larger concentrations shown on the 
right hand panel of Fig.[^ Subtracting the ionic contri¬ 
butions of Lii+m 

and N^ leads to better agreement 
with the expected linear behavior of X~^{T) over a wide 
temperature range (solid lines in Fig.[^. The summary 
of the obtained ^eff values and a comparison with the 
spin-only values are given in Tablejl] (extracted from the 
temperature range 150 K < T < 300 K). 

In the diluted systems /Teff signihcantly exceeds the 
spin-only value for the easy-axis systems with T = Fe and 
T = Ni in particular for field applied along the easy-axis. 
The diluted and the concentrated system of T = Fe show 
the largest enhancement of /igg when compared to the 
spin-only value and accordingly the magnetic anisotropy 




































TABLE I. Measured effective magnetic moments per tran¬ 
sition metal ion in Li 2 (Lii_a; ra;)N and the spin-only value 
calculated for in units of Bohr magneton. 
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T 

dilute 

H\\c H ±c 

concentrated 
H \\c 77 T c 

spin-only 

Mn 

5.5 

5.2 

3.6 

3.6 

4.9 

Fe 

6.7 

3.7 

6.5 

4.6 

3.8 

Co 

2.6 

3.1 

2.5 

2.7 

2.8 

Ni 

2.6 

2.0 

1.3 

1.3 

1.7 


observed in the isothermal magnetization measurements, 
M{H), is by far the largest among the different transi¬ 
tion metals (Fig.[^. In the easy-plane systems is only 
slightly enhanced. For the case of T = Mn this is directly 
reflected in the rather low magnetic anisotropy observed 
in M{H). Furthermore, the observation of smaller effec¬ 
tive moments for concentrated T = Mn when compared 
with the dilute system goes hand in hand with a further 
decrease of the anisotropy in M{H). In a similar fashion, 
the decrease of the effective moments for concentrated 
systems of both T = Co and T = Ni when compared 
to the dilute case is accompanied by a corresponding de¬ 
crease of the magnetic anisotropy in M{H). 

However, for dilute T = Co, a large magnetic 
anisotropy is apparent in M{H) (Fig.[^ even though the 
effective moments are only slightly enhanced when com¬ 
pared to the spin-only value. A more sophisticated anal¬ 
ysis is needed to describe the observed behavior. The 
input parameters for the employed analytical model are 
obtained from a simplified representation (Fig.[^) of the 
detailed electronic band structure. These were calculated 
earlier in the framework of local spin density approxima¬ 
tion (LSDA) for two different supercells which cor respond 
to Fe-concentrations of a; = 0.17^ and x = 0.50^^^®^, 
respectively. The calculated densities of states (DOS) 
projected on the 3d states of isolated Fe-atoms turn 
out to be fairly similar. In particular, a sharp peak 
in the minority spin channel that intersects the Fermi 
level appears in both cases. A further dilution of the 
Fe-concentration is not expected to cause tremendous 
changes of the projected Fe-DOS, since the Fe-atoms are 
already well separated. 

In our analytical model we consider only the DOS 
of the four 3d states with m = ±1 and m = ±2 and 
model them by Lorentzian-shaped peaks with a half¬ 
width w = 60meV (m: orbital quantum number). The 
TO = 0 states {3dz2) are not included in the model¬ 
ing since they are well below the Fermi level and have 
negligible effects on the MAE. The MAE depends on 
the crystal field splitting Ac (energy difference between 
TO = ±1 and TO = ±2 states), on the spin splitting As 
(energy difference between spin-up and spin-down), and 
on To facilitate the calculation of the MAE for 

T = Mn, Co, and Ni, only the Eermi level is shifted ac¬ 
cording to the number of 3d electrons leaving the band 
structure unchanged (rigid band approximation). Since 
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FIG. 3. (color online) a) Schematic density of states projected 
on the 3d states of isolated Fe atoms. The MAE in the analyt¬ 
ical model is determined by crystal field splitting (Ac), spin 
splitting (As), and peak width {w). For T = Mn, Co, and Ni 
the position of the Fermi level is shifted according to the 
number of 3d electrons (rigid band approximation), b) The 
magnetic anisotropy energy (MAE) per isolated T atom as a 
function of the number of 3d electrons calculated in the ana- 
l 3 rtical single-ion model is shown by the solid line. The MAE 
calculated by LSDA is shown for comparison (black squares). 
Inset are three single crystals on a millimeter grid and the lo¬ 
cal crystal structure of the transition metal with nearest and 
next-nearest neighbors. 


each degenerate state, and accordingly each Lorentzian in 
Eig.|^, can accommodate two electrons, the Fermi level 
either intersects the degenerate peaks or sits in the mid¬ 
dle, between two PDOS peaks (for an integer electron 
number). For the case of T = Mn^+, six 3d electrons 
have to be considered. Five of these occupy the ma¬ 
jority band (upper panel in Fig.j^) and the remaining 
one the to = 0 state of the minority band (lower panel 
in Fig.[^, the to = 0 state lies well below the Fermi 
level and is not shown). Accordingly, the Fermi level for 
T = Mni+ is located between the to = ±1 spin-up and 
the TO = ±2 spin-down states. For the case of T = Fe^+, 
one more electron has to be placed in the minority band. 
Therefore, one of the two to = ±2 spin-down states has 
to be occupied. Accordingly, the Fermi level intersects 
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the center of the m = ±2 spin-down peak. Adding one 
more electron shifts the Fermi level right between the 
m = ±2 and m = ±1 spin-down states corresponding 
to r = Co^’*' with eight 3d electrons. For T = Ni^+, 
one more 3d electron has to be considered and, similar 
to the case of T = Fe^"*", the Fermi level intersects the 
center of the m = ±1 spin-down peak. Furthermore, this 
allows for the calculation of the MAE as a continuous 
function of the band filling. It is worth mentioning that 
the spin orbit coupling is not included in the schematic 
band structure. 

The pairwise orbital susceptibility, and therefore the 
MAE, is now fully determined by w and the energy dif¬ 
ference between the Fermi level and the involved orbitals. 
The later one includes the effect of spin splitting and 
crystal electric field splitting. The pair susceptibility is 
proportional to the term w/[(eF — where ep 

is the Fermi energy, the DOS peak position with a 
= {spin-up, spin-down}!^. For T = Fe the m = ±2 
spin-down state is located right at the Fermi level, that 
is ep = e\. This leads to a resonance-like enhancement 
of the pair susceptibility that is only limited by the peak 
width w. The MAE for T = Fe is therefore dominated by 
the contribution of the m = ±2 spin-down state. This re¬ 
sult is in full analogy to the LSDA-based calculation, that 
revealed a splitting of this state cause d by sp in-orbit cou¬ 
pling which leads to large MAE values^i^ For T = Mn 
and Co no such resonance of the pair susceptibility oc¬ 
curs and the resulting MAE, which is determined by the 
sum over all orbital pairs, turns out to be negative. 

The MAEs calculated in our analytical model are 
shown in Fig.|^ as solid line (positive values correspond 
to an easy-axis, negative ones to an easy-plane system). 
Values obtained by LSDA methods!^ are given for com¬ 
parison (square data points). The MAE of Mn, Fe, Co, 
and Ni correspond to the respective integer number of 
3ci electrons. Most remarkably, this simplified model is 
sufficient to capture all basic features of the mag netic 
anisotropy even though the ’exact’ band structure h^ l ^^ l 
differs in several details from the schematic representa¬ 
tion shown in Fig.[^. The largest MAE is calculated for 
T = Fe. Given the simplicity of the model, the calculated 
value of 15 meV is in reasonable agreement with our ex¬ 
perimental result of 13meV for the dilute and 27meV 
for the more concentrated case (estimated from the lin¬ 
early extrapolated anisotropy field and the measured sat¬ 
uration magnetization) and with the LSDA result. The 
MAE calculated for T = Co is significantly smaller and 
of opposite sign. The smallest MAE is calculated for 
T = Mn. Both, sign and relative magnitude of the calcu¬ 
lated MAE fit well to the experimentally observed mag¬ 
netic anisotropy (Fig.[^. One exception is observed for 
the T = Ni system. The single-ion model gives a MAE of 
about half of the T = Fe value that results mainly from 
an |mp dependence of the MAE^. However, the exper¬ 
iment suggests that the MAE for concentrated T = Ni 
is reduced by more than an order of magnitude when 
compared to T = Fe. This discrepancy is mainly caused 


by an underestimation of the DOS peak width of the 
\m\ = 1 spin-down states that turns out to be larger 
than the |to| = 2 state (the MAE decreases with increas¬ 
ing band width, roughly following 1/w). Such a peak¬ 
broadening corresponds to an increasing delocalization 
which is indeed manifested in the decreasing electrical 
resistivity of Li 2 (Lii_a;Nia;)N for x > 0.8^^. Furthermore, 
spin-splitting and crysta l field splitting for T = Ni are 
smaller than for T = Fe h^ l ^^ l Adjusting the schematic 
band structure accordingly leads to a better agreement 
between single-ion model and experiment. There is also 
good agreement between our experimental results and 
recent calculations of the MAEs based on LSDA calcula¬ 
tions^. 

The orbital magnetic moment and the associated large 
MAE are likely direct consequences of the local symmetry 
of the transition metal. This actual linear, two-fold co¬ 
ordination between the nearest-neighbor nitrogen atoms 
gives rise to an effective, linear molecule (N-T-N). And as 
such, it is not subject to the Jahn-Teller effect^ which 
is driven by lifting the orbital degeneracy. Therefore, a 
quenching of the orbital magnetic moment by a lattice 
distortion does not take place in Li 2 (Lii_a; ra;)N. Within 
this symmetry the changing of the transition metal gives 
rise to a dramatic change of both experimental and cal¬ 
culated anisotropies. It remains to be seen whether this 
behavior is generic to linear complexes or restricted to 
the special case of Li 2 (Lii_a; T 2 ,)N. Further indications 
for the relevance of a linear arrangement to the formation 
of orbital magnetic moments and large MAE in 3d transi¬ 
tion metals can be found in seemingly unrelated systems: 
Ad-atoms on surfaces!^ (diatomic molecules built from 
substrate oxygen and adsorbed cobalt) and some linear 
transition metal complexes^ indeed show significant or¬ 
bital contributions to the magnetic moment. 

In summary, we found significant orbital contributions 
to the magnetic moment of the transition metals Mn, Fe, 
Co, and Ni substituted in LiaN. In accordance, large mag¬ 
netic anisotropies are observed. A sharp peak of the DOS 
which is intersected by the Fermi level gives rise to the 
uniaxial magnetic anisotropy of Fe and Ni. Even though 
this is not the case for Mn and Co, the latter one does also 
show a sizable magnetic anisotropy which is, however, of 
easy-plane type. This behavior can be described in an 
analytical, single-ion model based on only three parame¬ 
ters: crystal field splitting, spin splitting, and peak width 
of the DOS. Based on these considerations, it could be 
possible to identify, or even design, further magnetically 
ordered transition metal compounds with large orbital 
magnetic moments and magnetic anisotropy without re¬ 
lying on detailed band structure calculations and exces¬ 
sive computer power. 
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